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Abstract 

We deal with the problem of determining an inclusion within an elec- 
trical conductor from electrical boundary measurements. Under mild a 
priori assumptions we establish an optimal stability estimate. 

1 Introduction 

In this paper we deal with an inverse boundary value problem which is a special 
instance of the well-known Calderon's inverse conductivity problem 0. Given 
a bounded domain in M", n > 2, with reasonably smooth boundary, an open 
set D, compactly contained in Q, and a constant fc > 0, fc ^ 1, consider, for any 
/ G H^/^{dQ), the weak solution u G iJ^(ri) to the Dirichlet problem 



where xd denotes the characteristic function of the set D. We will denote by 
Ad ■■ H^^^idn) H-^/^{dn) the so called Dirichlet-to-Neumann map, that is 
the operator which maps the Dirichlet data onto the corresponding Neumann 
data f^iQfj- The inverse problem that we examine here is to determine D when 
A/3 is given. 

In '88 Isakov ^ proved the uniqueness, the purpose of the present paper 
is to prove a result of stability. In fact we prove that, under mild a priori 
assumptions on the regularity and on the topology of D, there is a continuous 
dependence of D (in the Hausdorff metric) from A^i with a modulus of continuity 
of logarithmic type, see Theorem l2.2l below. Let us stress that, indeed, this rate 
of continuity is the optimal one, as it was shown by examples in the recent paper 
|DC-R| by the second author and Luca Rondi. 

We wish to mention here a closely related, but different, problem which at- 
tracted a lot of attention starting from the papers of Friedman |0 and Friedman 
and Gustafsson [F-Cxj . That is the one of determining D when, instead of full 
knowledge of the Dirichlet-to-Neumann map, only one, or few, pairs of Dirich- 
let and Neumann data are available, see |A-I| . |I2] for extended bibliographical 
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accounts. Unfortunately, for such a problem, the uniqueness question, not to 
mention stability, remains a largely open issue. 

Let us illustrate briefly the main steps of our arguments. We must recall 
that Isakov's approach to uniqueness is essentially based on two arguments 

a) the Runge approximation theorem, 

b) the use of solutions with Green's function type singularities. 

Also here we shall use singular solutions, and indeed we shall need an accurate 
study of their asymptotic behavior when the singularity gets close to the set of 
discontinuity dD of the conductivity coefficient l + {k—l)xD in see Propo- 

sition lT^ On the other hand, it seems that Runge's theorem, which is typically 
based on nonconstructive arguments, (Lax, Kohn and Vogelius IK-Vp is not 
suited for stability estimates and therefore we introduced a different approach 
based on quantitative estimates of unique continuation, see Proposition 

In Section 2 we formulate our main hypotheses and state the stability result. 
Theorem 12. 21 In Section 3 we prove Theorem 12 . 21 on the basis of some auxiliary 
Propositions, whose proof is deferred to the following Section 4. 

2 The main result 

Let us introduce our regularity and topological assumptions on the conductor 
fl and on the unknown inclusion D. To this purpose we shall need the following 
definitions. In places, we shall denote a point x £ M" hy x = {x',Xn) where 
x' e Xn e R. 

Definition 2.1. Let Q be a bounded domain in M". Given a, < a < I, we 
shall say that a portion S of dfl is of class C^'" with constants r, L > if, for 
any P £ S , there exists a rigid transformation of coordinates under which we 
have P — and 

Q. n B-{0) ^ {x e B- : x„> (fix')}, 

where if is a C^'" function on 3^(0) C W"'^^ satisfying v(0) — |V</?(0)| — and 
ll<y5||ci.-(B-(o)) < Lr. 

Definition 2.2. We shall say that a portion S of dfl is of Lipschitz class with 
constants r, L > if for any P £ S, there exists a rigid transformation of 
coordinates under which we have P — and 

Q. n Br{0) = {x e Br : Xn> ^{x')}., 

where ip is a Lipschitz continuous function on -Br(O) C K"^"*^ satisfying Lp{0) = 
and |l.^|lco,i(B-(o)) < Lf. 

Remark 2.1. We have chosen to scale all norms in a such a way that they 
are dimensionally equivalent to their argument. For instance, for any ip G 
C^'"(B-(0)) we set 

ll'^llci.°(S-(0)) = I|¥'I|l~(B-(0)) +^l|V(^||L-(i3-(0)) +r'^^°'\^'P\a,B-{0)- 

For given numbers r, Af , (5, i > 0, < a < 1, we shall assume 
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(HI) the domain Q. satisfies the following conditions 

(2.1) \n\ < Mr"", 
where | • | denotes the Lebesgue measure of f2, 

(2.2) dO. is of class C^'" with constants f, L, 

(H2) the inclusion D satisfies the following conditions 

(2.3) \ D is connected, 



(2.4) dist(£>,aO) > (5, 



(2.5) dD is of class C^'" with constants r, L. 

In the sequel wc shall refer to numbers fc, n, r, M, (5, L, a as to the a priori 
data. We shall denote by Di and D2 two possible inclusions in both satisfying 
the properties mentioned. We shall denote by A/3., i = 1,2, the Dirichlet-to- 
Neumann map A/j when £) = £),. We can now state the main theorem. 

Theorem 2.2. Let n C M", n > 2, satisfy (HI). Let k > 0, k ^ I he given. 
Let Di and D2 be two inclusions in fl satisfying (H2). If, given e > 0, we have 

(2.6) II Adi - Ad^ I|£(hi/2,h-i/2) < e, 

then 

dnidDi,dD2) < w(e), 
where co is an increasing function on [0, +00), which satisfies 

uj{t) < C\ logij"'', for every 0<t<l 

and C,ri,C>0,0<r]<l, are constants only depending on the a priori data. 

Here d-n denotes the Hausdorff distance between bounded closed sets of M" 
and II • ||£(//i/2jj-i/2) denotes the operator norm on the space of bounded linear 
operators between H^/'^{dn) and H-'^/^{dn). 

Remark 2.3. It should be emphasized that in this statement the unknown in- 
clusion may be disconnected. 

Remark 2.4. Several variations of the above results could be devised with minor 
adaptations on the arguments. Just to mention one, an analogous result would 
be obtained if the Neumann-to-Dirichlet maps Nn^ are available instead of the 
Dirichlet-to-Neumann maps Ajj.. 
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3 Proof of Theorem [Zl 



Before proving Theorem 12. 21 we shall state some auxiliary Propositions, whose 
proofs are collected in the next Section 4. Here and in the sequel we shall de- 
note by Q the connected component of 51 \ {Di U -D2), whose boundary contains 
dfl, no = \g, = {x e Cn : dist(a;, Q) < r} and 827 = {x & M." : r < 
dist(a;,17) < 2r}. 

We introduce a variation of the Hausdorff distance which we call modified dis- 
tance. 

Definition 3.1. We shall call modified distance between Di and D2 the number 

(3.1) 1)2) = max < sup dist(a;,D2), sup dist(a;, Di) >. 

This notion is an adaptation of the one introduced in |A-B-R,-V] , which was also 
called modified distance. In order to distinguish such two notions, we call the 
present one, whereas the one in (A B-R-V^ was denoted by d„i. On the other 
hand, wc need to stress the common peculiarities: such modified distances do 
not satisfy the axioms of a metric and in general do not dominate the Hausdorff 
distance (see Section 3 in |A-B-R,-VI for related arguments). The following 
Proposition provides sufficient conditions under which dominates d-H. See 
|A-B-R-V) Proposition 3.6 for a related statement. 

Proposition 3.1. Let il be an open set in R" satisfying (HI). Let Di, D2 be 
two bounded open inclusions of VL satisfying (H2). Then 

(3.2) dn{dDi,dD2) < cd^{D^,D2), 
where c depends only on the a priori assumptions. 

With no loss of generality, we can assume that there exists a point O of dDi n 
drtn, where the maximum in the definition H3.1|l is attained, that is 

(3.3) d^=d^iDuD2)=dist{0,D2). 

As is well-known, the Dirichlet-to-Neumann map Ad associated to problem 
(Tmi . rr% is defined by: 

(3.4) < Adu, v>= [ (1 + {k - l)xD)yu ■ Vv, 

Jn 

for every u G H^{fl) solution to 1)1.1(1 and for every v e H^{il). Here < •, • > 
denotes the dual pairing between H-^^^{dn) and H^/'^{dVl). With a slight abuse 
of notation we shall write 

< 5>/ / gf da, 

J an 

for any / e H^/'^{d^) and g € H-^/^{dn). Let TD{x,y) be the fundamental 
solution for the operator div((l + {k — l)x£))V-), thus 

(3.5) div((l + (fc - 1)xd)VTd{; y)) = -S{- - y), 
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where y,w G R", 6 denotes the Dirac distribution . We shaU denote by Tu-^, 
Tr)2 such fundamental solutions when D = Di, D2 respectively. Recalling the 
well-known identity 

/ (l + (fc-l)xDi)VwrVu2- / (l+(fc-l)xDjViirVu2 = / mi [Aj^i -AdJu2, 
Jn Jn Jon 

which holds for every Ui e H^i^l), i = 1,2, solutions to when D = A 

respectively (see jI2j formula (5.0.4), Section 5.0), we have 

/,,(1 + (fc ~ l)xDjVrz,,(.,2;) • Vrz,, (.,«;) 

(3.6) - Ul + {k- l)xz3jVrz3, (•, 2/) • VTd2 i;w) 

= Ian^DA-,y)[-^Di - ■^D2](rD2{-,w))d(T, \/y,weCn. 
Let us define, for y,w E G 'J Cfl 

(3.7) SoAy^w) = (k-i) f vrz5,(-,zj)-vrz5,(.,u;), 

JDi 

(3.8) SD2{y,w) - (fc-1)/ VTDA;y)-^^D2{;w), 

JD2 

(3.9) f{y,w) = SdAv^w) - SD2{y,w). 
Thus (|3.6|l can be rewritten as 

(3.10) fiy,w)^ [ TDA-,yMD,-AD2WD2i-,w))da yy,weCTi. 

Jdn 

From now on we shall consider the dimension n > 3, since the case n — 2 
can be treated similarly through minor adaptations regarding the fundamental 
solutions. Up to a transformation of coordinates, we can assume that O, defined 
in H3.3|l . is the origin of the coordinate system. Let J^(0) be the outer unit 
normal vector to dflo in the origin O. Such a normal is indeed well-defined 
since we are assuming that O realizes the modified distance between Di and 
D2, therefore, in a small neighborhood of O, d^ln is made of a part of dDi, 
which is known to be C^'". We will rotate the coordinate system in such a way 
that i'{0) = (0, ... ,0, —1). Taking y = w ~ hv{0), with ft, > 0, we want to 
evaluate f{y, y) and Sdi [y, y) in term of h, for h small. Then, evaluating Sd2 in 
term of d^, we will get the stability estimate for the modified distance and thus, 
using Proposition 13. II for the Hausdorff distance. An important ingredient for 
evaluating / and Sdi is the behavior of the fundamental solution. We state 
now a proposition that collects all the results on Toi, i = 1,2, that we will 
need throughout the paper. For x = {x',Xn), where x' E M"^^ and Xn E M, 
we set X* = {x',—Xn)- We shall denote with x+ the characteristic function of 
the half-space {a;„ > 0} and with F^ the fundamental solution of the operator 
div((l + {k— l)x+)V-). If F is the standard fundamental solution of the Laplace 
operator, we have that (see for instance |A-FP| . Theorem 4) 

{ir(x, y) + j^^T{x, y*) for x„ > 0, y„ > 0, 

^T{x,y) forx„y„<0,, 

F(a;, y) - f^F(a;, y*) for x,, < 0, y,, < 0. 
The following Proposition holds. 
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Proposition 3.2. Let D C M" be an open set whose boundary is of class C^'", 
with constants r, L. 

(i) There exists a constant ci > depending on k, n, a and L only, such that 



(3.12) |V,ri3(x,y)| <ci|x-y|i-", 

for every x,y ^ M", 

(ii) There exist constants C2, C3 > depending on k, n, a and L only, such 
that 

(3.13) \TDix,y)-T+ix,y)\<^\x-y\' 

(3.14) \V,rDix,y)^V,T+{x,y)\ < - y\' 

for every x €z D D Br{0), and for every y — hv{0), with < r < Tq, 
0<h<ro, where Tq^ ( min {i(8L)-i/", i}) §. 

The next two Propositions give us quantitative estimates on / and Sdi when 
we move y towards O, along v{0). 

Proposition 3.3. Let il be an open set in K" satisfying (HI). Let Di,D2 be 
two inclusions in verifying (H2) and let y = hiy(0), with O defined in 
//, given e > 0, we have 

W^Di - A_D2||£(^i/2_^-i/2) < e, 
then for every h, < h < cf, where < c < 1, depends on L, 

Bh" 

(3.15) \f{y,y)\<C 



where Q < A < 1 and C,B,F > are constants that depend only on the a priori 
data. 

Proposition 3.4. Let be an open set in R" satisfying (HI ). Let Di, D2 be two 

inclusions in verifying (H2) and y = hv{0). Then for every h, < h < ro/2, 

(3.16) \SDdy,y)\ > ci/i'"" ~ C2dl'^^ + C3, 

where c\ , C2 and C3 are positive constants only depending on the a priori data. 
Here To is the number introduced in Provosition \H. M 

Now we have all the tools that we need to prove Theorem 12. 21 

Proof of Theorem\EE Let O G dDi satisfying (jSiSll, that is 

df,{Di,D2)^dist{0,D2)^d^. 

Then, for y ~ hv{0), with < /i < hi, where hi = min{d^, cr,ro/2}, using 
(TTlUl . we have 



(3.17) \s,Ay.y)\ < cj^^ _\^„_, _\^^_,dx . c^^^^\^2.-. 
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Using Proposition we have 

- \SD2iy,y)\ < ISdAv^v) - SD2{y,y)\ 



= \f{v,y)\<c-^. 

On the other hand, by Proposition 13 . 41 and (|3.17() 

\Sd, (y, v)\ ~ v)\ > c,h^^'' - - M'"'"- 

Thus we have 



-Bh' 



That is 



^Bh" 



(3.18) > C5/i2-"(l-e^''^/i^), 

where A^n — 2~A, A>0. Let h = h{e) where h{e) = m.m{\\ne\~^ ,d^}, 
for < e < £i, with Si such that exp(-B| Ineil^/^) = 1/2. If d^, < |lne|-2V 
the theorem follows using Proposition 13. II In the other case we have 

< e^U'^-l"'' < exp ( ~ B\ \ne\'/^). 

Then, for any £, < £ < ei, 

{d, - h{e)f-^^ > C6/l(£)2-", 

that is 

(3.19) d^<C7\\ne\-^^ 
where 5 — 1/{2F). When £ > £i, then 

I In £ I ^ 2F 

dp < diamfJ < diamfi —. 

I ln£i|~2T' 

Finally, using Proposition 13. II the theorem follows. □ 

4 Proofs of the auxiliary Propositions 

We premise the proof of Proposition 13 . II with one lemma. 

Lemma 4.1. Let il be an open set in M" satisfying (HI). Let D be a bounded 
open inclusion of fl satisfying (H2). Then for every P £ dD, there exists a 
continuous path 7 in Q \ D with one end-point in P and the other on dfl, such 
that for every 267 

(4.1) |z-P| < cdist(z,L)), 

where c is a positive constant depending on the a priori data only. 
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Proof. Using Lemma 5.2 of |A-B-R-V] . (which adapted arguments due to Lieber- 
man JI]), we approximate dist{-,dD) with a regularized distance d such that 
d e C^ifl \ D) U Ci'"(r2 \ D) and the foUowing facts hold 

dist{x,dD) 

70 < r-- < 71, 

d(x) 

|V(i(y)| > ci for every ?/ G ^2 s.t. dist(y,9D) > bf, 

\\d\\l,a < C2r, 

where 70, 71, 6, ci and C2 are positive constants only depending on L and a. 
We define for < h < of, with a depending on L and a only, 

= {a; e f7 \ U : d(x) > h}. 

Arguing as in Lemma 5.3 of |A-B-R-V] . Eh is connected with boundary of class 
and 

(4.2) cih < dist{x,dD) <C2h, yxedEhdn, 

where ci, C2 are positive constants depending on L and a only. Let us fix 
P e dD. Let J^(P) be the outer unit normal to dD in P (we recall that dD is 
C^'"). Since (|4.2I) . there exists a point P' € such that P' = hv{P), where 
/i is a positive constant ci/i < h < 02/1. We denote by PP' the segment whose 
end-points are P and P'. Since Eh is connected, there exists a continuous path 
7' C Eh with one end-point P' and the other on dH.. Since 7' C iJ/i we have 
that for every x € 7', dist(a::, 913) > ch, where c is a positive constant. We then 
define 7 = 7' U PP' and the lemma follows. □ 

Proof of Proposition Let us fix P e dDi . We distinguish the two following 
cases. 

i) P e dDi n dg, 

ii) PedDi\ dg. 
If case i) occurs then, 

dist{P,dD2) = dist(P,l32) < d^,. 

Let us consider case ii). Let 7 be the continuous path constructed in Lemma 
10 from P to dn. Since P ^ dg, there exists z ejH dD2 n ^f^u- 

dist(z,P)i)< sup {dist(2:,i:ii)} < d^(i:)i,D2)- 

Thus 

|z-P| < cd^{Di,D2), 
where c > is the constant appearing in H4.1(l On the other hand 

dist(P,aPi2) < \z-P\. 

So we obtain that, for every P G dDi 

dist(P,9D2) < cd^(Di,D2). 
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Similarly one can show that for every Q e dD2 

dist(Q,aDi) < cd^{Di,D2). 

Then we conclude 

dH{dDi,dD2) < cd^{Di,D2). 

□ 

Proof of Provosition I.V. M Let us prove (i). 

Let us consider the case x Cz D and y G dD. The cases in which x,y £ D or 
x,y € CD are trivial. Let = |x — ?/|. Let c be a positive number less than 
1+2 V" ' distinguish the following two cases: 

a) dist(x, dD) < ch, 

b) dist(x, dD) > ch. 

Let us consider the case a). Let P G dD be such that |P — a;| = dist{x,dD). 
For every r > 0, let Qr{P) be the cube centered at P, with sides of length 2r 
and parallel to the coordinates axes. We have that the ball Br{P) is inscribed 
into Qr{P). In particular x e Qch{P)- On the other hand 

\P-y\>\y-x\-\P~x\>h{i~c). 

Then, due to our choice of c, \P — y\ > {2ch)^pn, that is y ^ Q2ch{P)- Thus 
div,((l + (/c - 1)xd)V,Td{z. y)) =0 in Qa.JP) 

and for the piecewise C^'" regularity of F/j, proved in |DB-E-F| . see also |L- V| . 
we have 

(4.3) l|VFz,(-,2;)|U»(Q,,(p)) < |-||ri,(-, y)|U~(Q3,jp)), 

where ci depends on L, fc, n and a only. Using the pointwise bound of F^j with 
F (see |L-^{- W| V we have 

(4.4) ||Fc(-,y)|U^(Q3^jp)) <C2 (^yj . 

where C2 depends on n and fc only. Hence, by (|4.3|l and (|4.4|l . we get 

(4.5) \V^TD[x,y)\ < ||VF,,(.,y)|U=.(Q^,(p)) <C3/i'"" = C3|a;-2/|i-", 

where C3 depends on L, k, n and a only. 
If case b) occurs, then Qah_{x) C -D. Hence 

\V^rD{x,y)\ < ||VFD(-,y)|lL~(Q (x)) < ^l'^^^'' 
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where C4, C4 depend on i, fc, n and a only. 
Let us prove (iz). 

Let us fix ri = min {i(8i)-i/"r, §}. Recalling Definition EH we have that 

dD n B-(0) = {x e B-(0) : x„ = fix')}, 

where e C^^°'{«"~^) satisfying ^(0) = |V(p(0)| = 0. Let 6* G C°°(R) be such 
that < 6* < 1, eit) = 1, for |<| < 1, 0{t) = 0, for \t\ > 2 and |^| < 2. We 
consider the following change of variables £, — ^{x) defined by 

It can be verified that, with the given choice of ri, the following properties of 
<i> hold 

(4.6) <f>(Q2n(0)) = Q2n(0), 

(4.7) <i>(g,,(o)ni?) = Q+(o), 

(4.8) c-^i-xal < |$(a;i)-$(x2)| <c|xi-X2|, Va;i,a;2GM", 

(4.9) \^ix)-x\<^\x\^+°', VxeR", 

(4.10) \D<i>{x)~I\<^\x\'', VxgM", 

where (0) — {x G Qn (0) : x„ > 0} and c > 1 depends on L and a only. <i> 
is a C^'" diffeomorphism from R" into itself. Let us define the cylinder as 

Cr, ={a;eK" : \x'\ <ri,|a;„| <n}. 

For x,y G C^, we have that Td{£,,v) = ^Dix,y), where ^ = ^{x), rj — $(y), is 
solution of 

(4.11) dive((l + {k~ l)x+)i?(e)V4f c(e, V)) = -S{( - V), 

where B = with J = ||(<i>~^(0)- We observe that B is of class C" and 

i3(0) — I. Let us consider 

R{x,y) = rD(a;,y) - r+(x,y), 

where we keep the notation x, y to indicate ^, rj. By the properties of r+ and 
by H4.11|l . R satisfies 

div,((l + (fc - l)x+)V,i?(x, y)) = div,((l + (fc - l)x+)(/ - B)V,fD{x, y)). 

Let L > 0, depending on the a priori data only, be such that Vl C -82,(0). Thus 
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using the fundamental solution r+ we obtain 



-Rix,y)^ / (l + (fc-l)x+)(S-/)V,r+(2,zj).V,rc(z,x)dz 



+ (i + (fc-i)x+) 



dij{z) 



B^(o)nCri 



(1 + (A; - - /)V,r+(z,y) • V 



7. 



+ / (1 + (fc - 1)X+)(S - /)V3r+(z, y) • V 

/s^(0)xC,, 



9-Bi,(0) 



^(2;, ^)^^ (2^,2/) - r+(z,?;)^(a;,z) 



da{z) 



For |a;|,|?/|<ri/2, the last two integrals are bounded. Using (|3.12|l we obtain 



\R{x,v)\ < c[l+ \znx - z\'-^\y - z\'-^'dz 
Jc„ 



where c depends on L, a, k and n and 



h = I |zr|a;-z|i-"|2;-zp-"d^, 
I {\z\<4h}nc^^ 



Now 



l2= I |zr|a;-zr"|2;-zr"dz. 
I {\z\>Ah}nCr--^ 



h < 



'|-!«|<4 



|to|<4 



where /i = |x — ?/| and 



J\W\<4: 



and = x/h and ry = y/^- From standard bounds (see, for instance, [M] 
Chapter 2, Section 11), it is not difficult to see that 

Fi£., v) 5: const. < 00, 
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for all ^, e M", |$ - = 1. Thus 

/i <c|x-yr+2-". 

Let us consider now l2- Since \y\ = —yn < \x — y\ — h, we can deduce \z\ < 
^\y — z\ and \z\ < 2\x — z\ and thus obtain that 



J\z\>Ah 

Then we conclude 

(4.12) \R{x,y)\<c\x-yr+^-\ 

for every \y\ < ri/2, where c depends on L, a, k and n only. Let us go 
back to the original coordinates system. We observe that if x € ^^^{B^^^^{0)) 
and y = e„?/„, with y„ G (— ri/2,0) then |$(a;) — a;| is bounded by c\x — yl^"*"". 
Namely, since $(x) • y < and $(y) — y, by H4.8|l we have 

(4.13) c~^\x\ < |$(x)| < |$(x) - 2/1 < c\x - y\. 
On the other hand, by (|4.9(l and (|4.13() 

(4.14) |<i>(x) - x| < <^\x- 
We have 

-R(2^, y) =^D(x,y)- r+ (a;, y) 
= Td{x, y) - r+(x,y) + T+{^{x), ^{y)) - r+($(x), $(y)) 
= i?(<i>(a;), a>(y)) + r+(ci>(x), y) - T+{x, y). 

Using 113, ijlSl), and iBTHl we obtain 

\^D{x,y) - r+(a;,y)| 



< 


c 




yr 


h2-n 


< 


C 


lx- 


yr 


h2-ri 


< 


c" 

—a 


\x- 


yr 


l-2-n 



c' 



r 



I l+a J 1 — n 



where c" depends on fc, n, a and L only. We estimate now the first derivative 
of R. To estimate the first derivative of R let us consider a cube Q C B^^^^{x) 
of side cri/4, with < c < 1, such that x S dQ. The following interpolation 
inequality holds: 

l|Vi?(-,y)|U^(Q) < c||i?(-,y)||i-^(Q)|Vi?(.,y)|i,Q, 
where 5 = ^ depends on L only and 

I Vi?|a,Q - sup 1 _ ,|a ■ 



12 



Since, from the piecewise Holder continuity of VFu see H4.3|l . and also of Vr+, 
see (|3.11|) . we have that 

|vi?(-,y)U,Q < |Vrz3(.,y)U,Q + |vr+(.,y)U,Q < c/i-"+i-", 

where c depends on L only, thus wc conclude 

|V,i?(a;,zj)| < ^/j("+2-«)(i-^)/,(-"+i-n)5 = 

2 

where 77 = 7^ — . Thus 

(4.15) |V,^(x,2/)|<|^|x-y|"+i-", 

2 

where 77 — and c depends on L only. Concerning we have 
|V,r+($(a;),y)-V,r+(x,t/)| 
= |i?$(x)^Vr+(., y)|^(,) - V,r+(a;, y)| 

< |p$(x)^-/)vr+(.,y)|$(,)| 

+ |Vr+(-,y)|<i,(^) - Va;r+(a;,y)| 

< ^||Vr+(.,y)|U»(Q„^)|a: - ci>(.x)| + |Vr+(.,z/)U,Q|<i>(a;) - x\" 



where c depends on k, n, a and L only. □ 

Proof of Provosition \S. !A Let us fix y £ 827 and let us consider /(y, •)■ We have 
that 

(4.16) A^f{y,w)={) in CUd- 
For w € S2r, by (EiHl, (|^TI7|> and ifTTa we have 

(4.17) |/(y,u;)| <C(r,L,M)||Ac, -Az,,|| =e. 

Let us now estimate /(y, w) when w € Q. We define Cj'* = {a; € : dist(a;, ^l^i) > 
ft-}. For every w € , we, have that 

|5z?i(y,^«)| < \k-l\l |V,Fz5,(x,f)||V,Fc,(x,u;)|da: 

(4.18) ^ c /" |a;- < c/i^-". 

Similarly jS'^i^ (y, u;)| < ch^^". Then we conclude that 

\f{y,w)\<ch^-^ ing^. 
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At this stage we shall make use of the three spheres inequality for supremum 
norms of harmonic functions v, see for instance |K-M| . For every li, I2, 
1 < h < I2 and for every x S C/ U ^2- U fi- there exists t G (0, 1] , depending only 
on ^1, I2 and n such that 

We apply it for «(•) = f{y,-) in the ball Br(x), where x € S2r be such that 
dist(x,r) = r/2, where T = {a; e M" : dist(a;,n) = f } C 5^2^, = 3r = 3r/2 
and ^2 = 4r = 2r, then we obtain 

(4.19) \\f{y, OIlL-CSa^/sC^)) ^ 11/(1/' •)llL~(B-/2(3;))ll/(y' ') II ^"^"'(5,^(5)) ■ 

For every W G G'^, we denote with 7 a simple arc in U f2- U 52- joining x 
to W. Let us define {xi}, i = 1, . . . ,s as follows xi — x, Xi^i = 'yiU), where 
ti = max{t : \j{t) — Xi\ = r} if \xi — w\ > f, otherwise let i = s and stop the 
process. By construction, the balls BY/2{xi) are pairwise disjoint, \xi+i —Xi\ =r 
for i = 1, . . . , s — 1, |xs — Wl < r. For (|2.1|l . there exists /? such that s < (3. 
An iterated application of the three spheres inequality H4.19|l for f{y,-) (see for 
instance (A-B-R-V) pg.780, ^A-DB; Appendix E) gives that for any r, < r <f 



(4.20) ||/(y,-)llL^(s./.(^)) < ll/(y,-)IIL^(s,,,(^))ll/(y,-)lli^^(e)- 

We can now estimate the right hand side of H4.20|l by H4.17|l and (|4.18|) and 
obtain, for any r, < r < r 

(4.21) \\f{y, .)||l»(b.,,(^)) < c{h'-n'-^'e^' < c{h'-^)^s^ 
where /3 — and A = 1 — (3. Let O E dDi, as defined in H3.3|l . that is 

d{0,D2) = d^,{Di,D2). 

There exists a C^'" neighborhood [/ of O in dilo with constants r and L. Thus 
there exists a non-tangential vector field i?, defined on U such that the truncated 
cone 

(4.22) C{O,l^{O),0,r) = LeW : > cos 0, \x-0\<r 

[ \x-0\ 

satisfies 

CiO,^O),0,r)cg, 
where 9 = arctan(l/L). Let us define 



Ai = min 



h =^ arcsm 



r r 



1 + sin 6 3 sin t 



smf 



4 

wi = O + Xiv, 
pi — Ai sin6'i. 
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We have that Bp,{wi) C C{0,l){0),ei ,r), B jp.jwi) C C{0,li{0),9,r). Let 
w = wi, since pi < r/2, we can use (|4.21|) in the baU Bp^{w) and we can 
approach O G dDi by constructing a sequence of baUs contained in the cone 
C{0,v{0),ei,r). We define, for k>2 

^ . , , 1 — sin 6*1 

Wk=0 + Xkv, Xk^X^k^i, Pk = XPk-i, with X = :;— — 

1 + sm 01 

Hence pk = x'^^Vi^ = X^^^^i and 

Bp^^^{wk+i) C Bp.^^{wk) C Bp^^{wk) C C{0,v,e,r). 

Denoting d[k) = \wk ~ 0\ — pk ~ — pfc, we have d[k) = x'^~"^c^(l), with 
(i(l) = Ai(l — sinfl). For any r, < r < d{l), let A;(r) be the smahest integer 
such that d{k) < r, that is 

|logx| |logx| 

By an iterated apphcation of the three spheres inequahty over the chain of balls 

Bp,{wi), . . .,Bp^^^^^{wk(r)), we have 

(4.23) < c(/ii-")^e'3^'"^'", 

for < r < cr, where < c < 1 depends on L only. 



Let us consider now f{y, w) as a function of y. First observe that 

Ayf{y,w)^0 inCr^n, for all w G Cfi^). 
For y,w £ , y 7^ using 13.12|l . we have 

\SDAy,^)\<c I \x~y\'-^\x~w\'-''dx<ce-'\ 

Similarly for ■ Therefore 

\f{y, w)\ < ch"^-^"- with y,we . 
Finally, for y e 827 and w £ G^ , using 14.23|l . we have 

Proceeding as before, let us fix w G such that dist(zi;, Qfi^i) = h and y G 827 
such that dist(y, F) = r/2. Taking r = r/2, = 3r, I2 = 4r, j/i = O + \iv and 
using iteratively the three spheres inequality, we have 
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where t and s are the same number established previously. Therefore 

li~(B-/2(yi))) 



where 7 = r^, with as before, so < 7 < 1, and A' = At'' + 1 — 7. Once 
more, let us apply iteratively the three spheres inequality over a chain of balls 
contained in a cone with vertex in O and we obtain 

(4.24) < c(/,--)-'(-^^'^'-^)(s^^^'^'-^)-^""-\ 

Now, from H4.24(l . choosing y — w — hv{0), where ^{O) is the exterior unit 
normal to d^D in O, we obtain 

(4.25) |/(y,y)|<c/.-4"(£^^^'^'~^)^^^-'^'-\ 

where A" — —{2 — 2n)[3A' > 0. We observe that, for < /i < cr, where 
< c < 1 depends on L, k{h) < c\ logh\ ~ — clog/i, so we can write 

^k{h) _ g-cloghlog(T) _ ^-clogr _ ^c|logT| _ j^F ^ 

with F = c\ logrj. Therefore 

my,y)\ < h-^"e^^'''' 

^ ^-A" log h^Br "C'' log e 



^-A" log h+B'h'' loge 



Then in H4.25|l we obtain that 



□ 

Proof of Proposition \3.4\ Let us consider y ~ hi>{0), where v^O) is the exterior 
outer normal to dflo in O with O defined as in H3.3|l . < h < Tq, where fo is 
the number introduced in Proposition 13.21 and x G Di such that |a; — j/| < r, 
with < r < ro- Let us first observe that since O £ dDi and x G Di, for Tu-^ 
we have the asymptotic formula (|3.14() . which says that 



yx^Diix,y) - Va;r+(a:,y) 



<ci\x-y\'-"+'. 



Furthermore, since we are in the situation in which x G Di and y ^ Di, for 
(13.11(1 . r+{x,y) = 2/{k + l)T{x,y), where T{x,y) denotes the standard funda- 
mental solution of the Laplace operator. Let us consider now T]j2{x,y). With 
our choice of O, a; and y, we know that y ^ D2 but we do not have any infor- 
mation on X, that is we do not know in which side of the interface dD2 it is. 
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Thus we have to distinguish different situations. 

If a; e Br{0) n Di n D2, then we have the asymptotic formula (|3.11() for 
and from Lemma 3.1 of IXI the following formula holds 



(4.26) 



VJ^Di (a;, v) ■ VJ:d2 {x, y) >c\x-y 



Consider now the case x £ {Di \ D2) H Br(0). In this region let us consider 
a smaller ball Bp{0) centered in O with radius p where Q < p < d^. Since the 
definition of we have Bp C] D2 = ^- If x and y are in Bp{0), we have 



(4.27) 



A,{rDAx,y)-nx,y)) =0 inBp(O), 



1 [^D2{x,y)~T{x,y)]^gg^^^^ 
Thus by the maximum principle 

(4.28) |rz5,(a:,y)-r(x,y)| < cip2-" Vx,2/eBp(0), 
and by interior gradient bound 

(4.29) IV^ToA^.y) - V,T{x,y)\ < cap^"" Va; G 5^/2(0), Vy G Bp(0). 
Thus, using Lemma 3.1 of 0, in Bp/2{0) we obtain the formula 

(4.30) V^Td, (x, y) ■ V^Td, {x, y) > caja: - yp-^" _ 04^'-"". 

Let us consider /i < ro/2 and 5, (0) = {a; e M" : |a; - 0| < r}, with < r < tq- 
Then we have 

\SDAy,y)\ 

= |fc-l| j VTz^i • Vrz5,dx+ j VTz^i • VTz^.dx 

DinB^(O) I>i\B,.(0) 



> |fc-l| y VLi,, • VF^^dx - |/c-l| 

DinB,.(0) Di\B,.(0) 

The first term can be estimated as follows 

DinBr(o) 

\/Td, ■ VFD.dx + / VTd, ■ VFc.dx 

{DinD2)nBr{0) J{Di^D2)nBr{0) 

\/Td, ■ VFd.cJx + / VFu, • VFu^da; 

(I3inD2)nB,(0) J{Di^D2)nBp(0) 



> 



[(Z5i^I52)nS,(0)]^Bp(0) 



VFi3, • VFi3,dx 
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In conclusion, choosing p = and using (|4.26f) and (|3.12l) we obtain 

\SDiy,y)\ > ci I \x-y\'-'^dx 

-C2 J |x-yr-"|x-y|i-"dx 

-C3 J \x^y\'--\x-y\'--dx 

Di^B^O) 
> C4h^-"- - Csd^"^" - C7. 

□ 
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